Stress Knowledge Map: A knowledge graph resource for systems biology analysis of plant stress responses.

Stress Knowledge Map (SKM, https://skm.nib.si) is a publicly available resource containing two complementary knowledge graphs describing current knowledge of biochemical, signalling, and regulatory molecular interactions in plants: a highly curated model of plant stress signalling (PSS, 543 reactions) and a large comprehensive knowledge network (CKN, 488,390 interactions). Both were constructed by domain experts through systematic curation of diverse literature and database resources. SKM provides a single entrypoint for plant stress response investigations and the related growth trade-offs as well as interactive exploration of current knowledge. PSS is also formulated as qualitative and quantitative models for systems biology, and thus represents a starting point of a plant digital twin. Here, we describe the features of SKM and show, through two case studies, how it can be used for complex analyses, including systematic hypothesis generation, design of validation experiments, or to gain new insights into experimental observations in plant biology.

Transcriptome-informed identification and characterization of Planococcus citri cis- and trans-isoprenyl diphosphate synthase genes.

Insect physiology and reproduction depend on several terpenoid compounds, whose biosynthesis is mainly unknown. One enigmatic group of insect monoterpenoids are mealybug sex pheromones, presumably resulting from the irregular coupling activity of unidentified isoprenyl diphosphate synthases (IDSs). Here, we performed a comprehensive search for IDS coding sequences of the pest mealybug Planococcus citri. We queried the available genomic and newly generated short- and long-read P. citri transcriptomic data and identified 18 putative IDS genes, whose phylogenetic analysis indicates several gene family expansion events. In vitro testing confirmed regular short-chain coupling activity with five gene products. With the candidate with highest IDS activity, we also detected low amounts of irregular coupling products, and determined amino acid residues important for chain-length preference and irregular coupling activity. This work therefore provides an important foundation for deciphering terpenoid biosynthesis in mealybugs, including the sex pheromone biosynthesis in P. citri.

Cell-type proteomic and metabolomic resolution of early and late grain filling stages of wheat endosperm.

The nutritional value of wheat grains, particularly their protein and metabolite composition, is a result of the grain-filling process, especially in the endosperm. Here, we employ laser microdissection (LMD) combined with shotgun proteomics and metabolomics to generate a cell type-specific proteome and metabolome inventory of developing wheat endosperm at the early (15 DAA) and late (26 DAA) grain-filling stages. We identified 1803 proteins and 41 metabolites from four different cell types (aleurone (AL), sub-aleurone (SA), starchy endosperm (SE) and endosperm transfer cells (ETCs). Differentially expressed proteins were detected, 67 in the AL, 31 in the SA, 27 in the SE and 50 in the ETCs between these two-time points. Cell-type accumulation of specific SUT and GLUT transporters, sucrose converting and starch biosynthesis enzymes correlate well with the respective sugar metabolites, suggesting sugar upload and starch accumulation via nucellar projection and ETC at 15 DAA in contrast to the later stage at 26 DAA. Changes in various protein levels between AL, SA and ETC support this metabolic switch from 15 to 26 DAA. The distinct spatial and temporal abundances of proteins and metabolites revealed a contrasting activity of nitrogen assimilation pathways, e.g. for GOGAT, GDH and glutamic acid, in the different cell types from 15 to 26 DAA, which can be correlated with specific protein accumulation in the endosperm. The integration of cell-type specific proteome and metabolome data revealed a complex metabolic interplay of the different cell types and a functional switch during grain development and grain-filling processes.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato.

BACKGROUND: Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV). RESULTS: We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes. CONCLUSION: We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Studying Cell Death Initiation Using a Digital Microscope.

Hypersensitive response (HR)-conferred resistance is an effective defense response that can be determined by the N resistance genes. HR is manifested as the formation of cell death zones on inoculated leaves. Here, a protocol for studying the rate of cell death initiation by imaging inoculated leaves in the time between the cell death initiation and the cell death appearance using a digital microscope is presented. The digital microscope enables a continuous imaging process in desired intervals, which allows an accurate determination of cell death initiation rate up to minutes exactly, as opposed to hours in traditional methods. Imaging with the digital microscope is also independent of light and can therefore be used during day and night without disturbing the circadian rhythm of the plant. Different pathosystems resulting in programmed cell death development could be studied using this protocol with minor modifications. Overall, the protocol thus allows simple, accurate, and inexpensive identification of cell death initiation rate.

Candidate pathogenicity factor/effector proteins of 'Candidatus Phytoplasma solani' modulate plant carbohydrate metabolism, accelerate the ascorbate-glutathione cycle, and induce autophagosomes.

The pathogenicity of intracellular plant pathogenic bacteria is associated with the action of pathogenicity factors/effectors, but their physiological roles for most phytoplasma species, including 'Candidiatus Phytoplasma solani' are unknown. Six putative pathogenicity factors/effectors from six different strains of 'Ca. P. solani' were selected by bioinformatic analysis. The way in which they manipulate the host cellular machinery was elucidated by analyzing Nicotiana benthamiana leaves after Agrobacterium-mediated transient transformation with the pathogenicity factor/effector constructs using confocal microscopy, pull-down, and co-immunoprecipitation, and enzyme assays. Candidate pathogenicity factors/effectors were shown to modulate plant carbohydrate metabolism and the ascorbate-glutathione cycle and to induce autophagosomes. PoStoSP06, PoStoSP13, and PoStoSP28 were localized in the nucleus and cytosol. The most active effector in the processes studied was PoStoSP06. PoStoSP18 was associated with an increase in phosphoglucomutase activity, whereas PoStoSP28, previously annotated as an antigenic membrane protein StAMP, specifically interacted with phosphoglucomutase. PoStoSP04 induced only the ascorbate-glutathione cycle along with other pathogenicity factors/effectors. Candidate pathogenicity factors/effectors were involved in reprogramming host carbohydrate metabolism in favor of phytoplasma own growth and infection. They were specifically associated with three distinct metabolic pathways leading to fructose-6-phosphate as an input substrate for glycolysis. The possible significance of autophagosome induction by PoStoSP28 is discussed.

Chloroplast redox state changes mark cell-to-cell signaling in the hypersensitive response.

Hypersensitive response (HR)-conferred resistance is associated with induction of programmed cell death and pathogen spread restriction in its proximity. The exact role of chloroplastic reactive oxygen species and its link with salicylic acid (SA) signaling in HR remain unexplained. To unravel this, we performed a detailed spatiotemporal analysis of chloroplast redox response in palisade mesophyll and upper epidermis to potato virus Y (PVY) infection in a resistant potato genotype and its transgenic counterpart with impaired SA accumulation and compromised resistance. Besides the cells close to the cell death zone, we detected individual cells with oxidized chloroplasts further from the cell death zone. These are rare in SA-deficient plants, suggesting their role in signaling for resistance. We confirmed that chloroplast redox changes play important roles in signaling for resistance, as blocking chloroplast redox changes affected spatial responses at the transcriptional level. Through spatiotemporal study of stromule induction after PVY infection, we show that stromules are induced by cell death and also as a response to PVY multiplication at the front of infection. Overall induction of stromules is attenuated in SA-deficient plants.

A mini-TGA protein modulates gene expression through heterogeneous association with transcription factors.

TGA (TGACG-binding) transcription factors, which bind their target DNA through a conserved basic region leucine zipper (bZIP) domain, are vital regulators of gene expression in salicylic acid (SA)-mediated plant immunity. Here, we investigated the role of StTGA2.1, a potato (Solanum tuberosum) TGA lacking the full bZIP, which we named a mini-TGA. Such truncated proteins have been widely assigned as loss-of-function mutants. We, however, confirmed that StTGA2.1 overexpression compensates for SA-deficiency, indicating a distinct mechanism of action compared with model plant species. To understand the underlying mechanisms, we showed that StTGA2.1 can physically interact with StTGA2.2 and StTGA2.3, while its interaction with DNA was not detected. We investigated the changes in transcriptional regulation due to StTGA2.1 overexpression, identifying direct and indirect target genes. Using in planta transactivation assays, we confirmed that StTGA2.1 interacts with StTGA2.3 to activate StPRX07, a member of class III peroxidases (StPRX), which are known to play role in immune response. Finally, via structural modeling and molecular dynamics simulations, we hypothesized that the compact molecular architecture of StTGA2.1 distorts DNA conformation upon heterodimer binding to enable transcriptional activation. This study demonstrates how protein truncation can lead to distinct functions and that such events should be studied carefully in other protein families.

pISA-tree - a data management framework for life science research projects using a standardised directory tree.

We developed pISA-tree, a straightforward and flexible data management solution for organisation of life science project-associated research data and metadata. pISA-tree was initiated by end-user requirements thus its strong points are practicality and low maintenance cost. It enables on-the-fly creation of enriched directory tree structure (project/Investigation/Study/Assay) based on the ISA model, in a standardised manner via consecutive batch files. Templates-based metadata is generated in parallel at each level enabling guided submission of experiment metadata. pISA-tree is complemented by two R packages, pisar and seekr. pisar facilitates integration of pISA-tree datasets into bioinformatic pipelines and generation of ISA-Tab exports. seekr enables synchronisation with the FAIRDOMHub repository. Applicability of pISA-tree was demonstrated in several national and international multi-partner projects. The system thus supports findable, accessible, interoperable and reusable (FAIR) research and is in accordance with the Open Science initiative. Source code and documentation of pISA-tree are available at https://github.com/NIB-SI/pISA-tree .

Transcriptional deregulation of stress-growth balance in Nicotiana benthamiana biofactories producing insect sex pheromones.

Plant biofactories are a promising platform for sustainable production of high-value compounds, among which are insect sex pheromones, a green alternative to conventional insecticides in agriculture. Recently, we have constructed transgenic Nicotiana benthamiana plants ("Sexy Plants", SxP) that successfully produce a blend of moth (Lepidoptera) sex pheromone compounds (Z)-11-hexadecen-1-ol and (Z)-11-hexadecenyl acetate. However, efficient biosynthesis of sex pheromones resulted in growth and developmental penalty, diminishing the potential for commercial use of SxP in biomanufacturing. To gain insight into the underlying molecular responses, we analysed the whole-genome transcriptome and evaluated it in relation to growth and pheromone production in low- and high-producing transgenic plants of v1.0 and v1.2 SxP lines. In our study, high-producing SxPv1.2 plants accumulated the highest amounts of pheromones but still maintained better growth compared to v1.0 high producers. For an in-depth biological interpretation of the transcriptomic data, we have prepared a comprehensive functional N. benthamiana genome annotation as well as gene translations to Arabidopsis thaliana, enabling functional information transfer by using Arabidopsis knowledge networks. Differential gene expression analysis, contrasting pheromone producers to wild-type plants, revealed that while only a few genes were differentially regulated in low-producing plants, high-producing plants exhibited vast transcriptional reprogramming. They showed signs of stress-like response, manifested as downregulation of photosynthesis-related genes and significant differences in expression of hormonal signalling and secondary metabolism-related genes, the latter presumably leading to previously reported volatilome changes. Further network analyses confirmed stress-like response with activation of jasmonic acid and downregulation of gibberellic acid signalling, illuminating the possibility that the observed growth penalty was not solely a consequence of a higher metabolic burden imposed upon constitutive expression of a heterologous biosynthetic pathway, but rather the result of signalling pathway perturbation. Our work presents an example of comprehensive transcriptomic analyses of disadvantageous stress signalling in N. benthamiana biofactory that could be applied to other bioproduction systems.

Corrigendum: Physiological and transcriptional responses to saline irrigation of young 'Tempranillo' vines grafted onto different rootstocks.

[This corrects the article DOI: 10.3389/fpls.2022.866053.].

Organelles and phytohormones: a network of interactions in plant stress responses.

Phytohormones are major signaling components that contribute to nearly all aspects of plant life. They constitute an interconnected communication network to fine-tune growth and development in response to the ever-changing environment. To this end, they have to coordinate with other signaling components, such as reactive oxygen species and calcium signals. On the one hand, the two endosymbiotic organelles, plastids and mitochondria, control various aspects of phytohormone signaling and harbor important steps of hormone precursor biosynthesis. On the other hand, phytohormones have feedback actions on organellar functions. In addition, organelles and phytohormones often act in parallel in a coordinated matter to regulate cellular functions. Therefore, linking organelle functions with increasing knowledge of phytohormone biosynthesis, perception, and signaling will reveal new aspects of plant stress tolerance. In this review, we highlight recent work on organelle-phytohormone interactions focusing on the major stress-related hormones abscisic acid, jasmonates, salicylic acid, and ethylene.

CRISPR/Cas9-mediated fine-tuning of miRNA expression in tetraploid potato.

MicroRNAs (miRNAs) are small noncoding RNAs, which modulate the abundance and spatiotemporal accumulation of target mRNAs at transcriptional and post-transcriptional levels and through that play important roles in several biological processes in plants. Here we show that in polyploid species, CRISPR/Cas9 system can be used for fine-tuning of miRNA expression, which can have broader range of applications compared to knock-out mutants. We established the complete pipeline for CRISPR-Cas9-mediated modulation of miRNA expression in potato. It consists of (1) design and assembly of dual sgRNA CRISPR/Cas9 constructs, (2) transient transfection of protoplasts following fast and efficient screening by high resolution melting analysis to select functional sgRNAs, and (3) stable transformation of potato explants with functional sgRNAs and selection of regenerated transgenic lines with desired mutations and desired miRNA abundance based on sequencing and RT-qPCR. We show that miRNA-editing using dual sgRNA approach results in different types of mutations among transgenic lines but also in different alleles of the same plant, which are target site-dependent. The most frequent were short deletions, but we also detected 1-nt insertions (T or G), deletions between two sgRNAs and larger deletions. miRNA abundance correlates with the frequency and type of introduced mutations, as more extensive mutations in more alleles result in lower miRNA abundance. Interestingly, some mutated loci can generate alternative miRNAs, now novel targets were however predicted for those. In all transgenic lines with Cas9 expression, we detected mutations, suggesting high efficiency of Cas9-editing. We confirmed the miRNA-editing efficiency of our optimised approach in two different potato genotypes and three different loci.

Toward learning the principles of plant gene regulation.

Advanced machine learning (ML) algorithms produce highly accurate models of gene expression, uncovering novel regulatory features in nucleotide sequences involving multiple cis-regulatory regions across whole genes and structural properties. These broaden our understanding of gene regulation and point to new principles to test and adopt in the field of plant science.

TGA transcription factors-Structural characteristics as basis for functional variability.

TGA transcription factors are essential regulators of various cellular processes, their activity connected to different hormonal pathways, interacting proteins and regulatory elements. Belonging to the basic region leucine zipper (bZIP) family, TGAs operate by binding to their target DNA sequence as dimers through a conserved bZIP domain. Despite sharing the core DNA-binding sequence, the TGA paralogues exert somewhat different DNA-binding preferences. Sequence variability of their N- and C-terminal protein parts indicates their importance in defining TGA functional specificity through interactions with diverse proteins, affecting their DNA-binding properties. In this review, we provide a short and concise summary on plant TGA transcription factors from a structural point of view, including the relation of their structural characteristics to their functional roles in transcription regulation.

Physiological and Transcriptional Responses to Saline Irrigation of Young 'Tempranillo' Vines Grafted Onto Different Rootstocks.

The use of more salt stress-tolerant vine rootstocks can be a sustainable strategy for adapting traditional grapevine cultivars to future conditions. However, how the new M1 and M4 rootstocks perform against salinity compared to conventional ones, such as the 1103-Paulsen, had not been previously assessed under real field conditions. Therefore, a field trial was carried out in a young 'Tempranillo' (Vitis vinifera L.) vineyard grafted onto all three rootstocks under a semi-arid and hot-summer Mediterranean climate. The vines were irrigated with two kinds of water: a non-saline Control with EC of 0.8 dS m-1 and a Saline treatment with 3.5 dS m-1. Then, various physiological parameters were assessed in the scion, and, additionally, gene expression was studied by high throughput sequencing in leaf and berry tissues. Plant water relations evidenced the osmotic effect of water quality, but not that of the rootstock. Accordingly, leaf-level gas exchange rates were also reduced in all three rootstocks, with M1 inducing significantly lower net photosynthesis rates than 1103-Paulsen. Nevertheless, the expression of groups of genes involved in photosynthesis and amino acid metabolism pathways were not significantly and differentially expressed. The irrigation with saline water significantly increased leaf chloride contents in the scion onto the M-rootstocks, but not onto the 1103P. The limitation for leaf Cl- and Na+ accumulation on the scion was conferred by rootstock. Few processes were differentially regulated in the scion in response to the saline treatment, mainly, in the groups of genes involved in the flavonoids and phenylpropanoids metabolic pathways. However, these transcriptomic effects were not fully reflected in grape phenolic ripeness, with M4 being the only one that did not cause reductions in these compounds in response to salinity, and 1103-Paulsen having the highest overall concentrations. These results suggest that all three rootstocks confer short-term salinity tolerance to the scion. The lower transcriptomic changes and the lower accumulation of potentially phytotoxic ions in the scion grafted onto 1103-Paulsen compared to M-rootstocks point to the former being able to maintain this physiological response in the longer term. Further agronomic trials should be conducted to confirm these effects on vine physiology and transcriptomics in mature vineyards.

Analysis of Virus Spread Around the Cell Death Zone at Spatiotemporal Resolution Using Confocal Microscopy.

The role of programmed cell death (PCD) in hypersensitive response (HR)-conferred resistance depends on the type of host-pathogen interaction and therefore has to be studied for each individual pathosystem. Here we present and explain the protocol for studying the role of PCD in HR-conferred resistance in potato plants in the interaction with the viral pathogen. As an experimental system, we use genotype Rywal, where the virus spread is restricted and HR PCD develops 3 days post potato virus Y (PVY) inoculation. As a control of virus multiplication and spread, we include its transgenic counterpart impaired in salicylic acid (SA) accumulation (NahG-Rywal), in which the HR-PCD occurs but the spread of the virus is not restricted. To follow the occurrence of virus-infected cells and/or virus multiplication outside the cell death zone, we use GFP-tagged PVY (PVY-N605(123)-GFP) which can be monitored by confocal microscopy. Any other plant-pathogen system which results in PCD development could be studied using a modified version of this protocol.

Fermentative Indole Production via Bacterial Tryptophan Synthase Alpha Subunit and Plant Indole-3-Glycerol Phosphate Lyase Enzymes.

Indole is produced in nature by diverse organisms and exhibits a characteristic odor described as animal, fecal, and floral. In addition, it contributes to the flavor in foods, and it is applied in the fragrance and flavor industry. In nature, indole is synthesized either from tryptophan by bacterial tryptophanases (TNAs) or from indole-3-glycerol phosphate (IGP) by plant indole-3-glycerol phosphate lyases (IGLs). While it is widely accepted that the tryptophan synthase α-subunit (TSA) has intrinsically low IGL activity in the absence of the tryptophan synthase ß-subunit, in this study, we show that Corynebacterium glutamicum TSA functions as a bona fide IGL and can support fermentative indole production in strains providing IGP. By bioprospecting additional bacterial TSAs and plant IGLs that function as bona fide IGLs were identified. Capturing indole in an overlay enabled indole production to titers of about 0.7 g L-1 in fermentations using C. glutamicum strains expressing either the endogenous TSA gene or the IGL gene from wheat.

PaintOmics 4: new tools for the integrative analysis of multi-omics datasets supported by multiple pathway databases.

PaintOmics is a web server for the integrative analysis and visualisation of multi-omics datasets using biological pathway maps. PaintOmics 4 has several notable updates that improve and extend analyses. Three pathway databases are now supported: KEGG, Reactome and MapMan, providing more comprehensive pathway knowledge for animals and plants. New metabolite analysis methods fill gaps in traditional pathway-based enrichment methods. The metabolite hub analysis selects compounds with a high number of significant genes in their neighbouring network, suggesting regulation by gene expression changes. The metabolite class activity analysis tests the hypothesis that a metabolic class has a higher-than-expected proportion of significant elements, indicating that these compounds are regulated in the experiment. Finally, PaintOmics 4 includes a regulatory omics module to analyse the contribution of trans-regulatory layers (microRNA and transcription factors, RNA-binding proteins) to regulate pathways. We show the performance of PaintOmics 4 on both mouse and plant data to highlight how these new analysis features provide novel insights into regulatory biology. PaintOmics 4 is available at https://paintomics.org/.

Plant X-tender Toolbox for the Assembly and Expression of Multiple Transcriptional Units in Plants.

Plant synthetic biology requires the design of plant expression vectors with multiple transcriptional units, which can be challenging. Here we describe the use of Plant X-tender toolbox complemented with a plant grammar implemented in GenoCAD for design, cloning and delivery of several transcriptional units into the plant genome. Plant X-tender toolbox consists of a set of plant expression vectors and the protocols for the most efficient cloning of multiple transcriptional units into this novel vector set. Together with the plant grammar implemented in GenoCAD, the presented strategy allows the users to quickly design genetic modules and assemble them into Plant X-tender expression vectors for in planta functional studies or synthetic biology applications.